APS/123-QED 



The Turbulent Magnetic Diffusivity Tensor for Time-Dependent Mean Fields 
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We re-examine the nature of the turbulent magnetic diffusivity tensor of mean field electrodynam- 
ics and show that an inconsistency arises if it is calculated via consideration of time-independent 
magnetic fields. Specifically, the predicted growth rate of the mean magnetic field is, in general, 
incorrect. We describe how the traditional expansion procedure for the mean electromotive force 
should be extended, and illustrate the consistency of this approach by means of a perturbation 
analysis for a mean magnetic field varying on long spatial scales. Finally we examine the magnitude 
of this new contribution to the magnetic diffusion for a particular flow. 
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Large-scale magnetic fields are observed in a vast range 
of cosmical bodies, from planets to stars and accretion 
discs. Understanding their generation by turbulent dy- 
namo action remains one of the most challenging prob- 
lems in astrophysics. Cosmical dynamos are often stud- 
ied within the framework of mean field electrodynamics, 
a tremendously elegant theory that describes the evolu- 
tion of a mean (large-scale) magnetic field in terms of 
transport coefficients determined from averaged small- 
scale properties of the flow and field. The main under- 
lying assumption of mean field electrodynamics is one 
of scale separation for the velocity and magnetic fields, 
which is nearly always understood as spatial scale sepa- 
ration. In this case, the velocity and magnetic field can 
formally be decomposed into mean and fluctuating parts, 



Averaging equation ([2]) over some intermediate scale 
a, satisfying l«n<L, leads to the following equations 
for the mean and fluctuating magnetic fields, 
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Vx{U xB ) + Vx£ + 7iV 2 B , (3) 



db 

— = Vx([/ x b) + Vx(u x B )+VxG + r]V 2 b; (4) 
at 

here ( ) denotes a spatial average, £ = (u x b) is the mean 
electromotive force (emf ) and G = (it x b) — (u x b) . 
Equation ([4]) can be written formally as 



C(b) = Vx(tixB ), 
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where the fluctuating quantities u and b vary on the 
(small) scale I, which may be regarded as a typical tur- 
bulent eddy size; the mean fields Uq and Bo vary on 
some much larger scale L 3> i. Mean field electrody- 
namics is, at heart, a kinematic theory; here the velocity 
U(x,t) is assumed to be known, with the magnetic field 
determined solely by the induction equation 



— = V x (U x B) + r 1 V 2 B, 
at 



(2) 



where rj is the magnetic diffusivity (here assumed uni- 
form). In this paper, in which we shall address a funda- 
mental issue in the interpretation of mean field electro- 
dynamics, we shall also adopt this kinematic approach. 
We shall not consider the dynamical back-reaction of the 
magnetic field on the flow via the Lorentz force, impor- 
tant though it is in natural dynamos. 



'Permanent address: Department of Applied Mathematics, Uni- 
versity of Leeds, Leeds LS2 9JT, UK; Electronic address: 
|d. w.h ughes@lee ds.ac.uk] 

T Permanent address: Centre for Mathematical Sciences, Univer- 
sity of Cambridge, Wilberforce Road, Cambridge CB3 OWA, UK; 
Electronic address: mrep@cam.ac.uk 



where £ is a linear operator. The traditional approach to 
tackling © [H, i] assumes that if there is no mean field 
Bq then there are no non-decaying solutions for 6; in 
other words, all solutions to the homogeneous equation 
C(b) = decay. In reality, at high magnetic Reynolds 
number Rm this is unlikely to be the case; in other words, 
small-scale dynamo action will ensue. This can lead to 
severe difficulties in the consistency of the mean field ap- 
proach, as discussed in Q. However, in this Letter, we 
shall not explore this particular issue. Instead, we shall 
assume, as is certainly the case for low Rm, that the 
fluctuating magnetic field b is indeed driven entirely by a 
non-zero right hand side of equation It then follows 
that b, and hence £, is linearly and homogeneously re- 
lated to the mean field Bo- Accordingly, one may posit 
an expansion for £ in terms of Bq and its derivatives. 
This is always written as 



Si = 



dR 



Oj 



dx k 



(6) 



where it is anticipated that the large spatial scale of Bq 
will lead to rapid convergence. 

The aim of mean field electrodynamics is first to cal- 
culate the coefficients of the tensors otij and (3ijk, and 
then to substitute for £ into ([3]) to obtain an equation 
determining the temporal evolution of a large-scale field. 
It is the inconsistency between these two steps that we 
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shall explore in this Letter. In general, expression ([6]) 
should also contain temporal derivatives of B , although, 
as pointed out in [l[, these may be replaced by spatial 
derivatives by using back-substitution from (J3]). How- 
ever, and this is the crux of our argument, although this 
manipulation is formally correct, it leads to an inconsis- 
tency in the determination of the evolution of the mean 
field when the expression for £ is substituted from Q 
into fl3j). 

The difficulty arises as a consequence of the natural 
means of determining the tensors and /3j,-/- in expres- 
sion ([5]). Given a turbulent flow u(x,t), the components 
ctij can be determined by evaluation of the mean emf 
£ following the imposition of three independent, steady, 
spatially uniform fields Bq. Having calculated all the 
components of the onj tensor, the components /Jy^ can 
subsequently be determined by consideration of the emf 
following the imposition of independent, steady fields 
having a uniform gradient. Substitution into ([3]) then 
leads to an expression for the growth rate of the mean 
field in terms of the a and (3 tensors. However, in general, 
this expression will not provide the correct description of 
the evolution of a large-scale field. The problem arises 
because it is inconsistent to calculate fiijk from a time- 
independent spatially-dependent field, since, in reality, a 
field with spatial dependence will vary with time. To 
obtain a correct representation for the growth rate, pre- 
mature back-substitution for the temporal derivatives of 
the mean field must be avoided. Instead of ©, the fol- 
lowing expansion should be used for the mean emf: 



£j = ctijBoj + Tij — — — h pijk - 



(7) 
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where the coefficients Qy and [3ijk are identical to those 
in expression ([B]) and can be calculated precisely as de- 
scribed above, using time-independent magnetic fields. 
The new tensor is to be determined via evaluation of 
£ after imposition of a spatially uniform magnetic field 
that increases linearly with time. Substitution of the ex- 
pression for £ from ([7]) into ([3]) (ignoring both the mean 
flow and the molecular diffusion term, neither of which 
is important for the argument we are advancing here) 
yields an equation of the form 



dB 0l 
dt 
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Making use of the fact that the expression for the emf is 
assumed to be a rapidly convergent series, we may now 
back-substitute for the time derivative of Bq on the right 
hand side of ((U using just the leading order terms, thus 
yielding 

dB Ql _0_ 



Of 



dxj 

r —( r u r dB ° r< 
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The new term involving Tk m has one spatial derivative 
and hence is of the same order as that involving ftkmn- 



For simplicity, if we consider the case when the com- 
ponents of the a, (5 and T tensors are constants (more 
generally, they could be functions of the slow spatial vari- 
ation), it can be seen that the coefficient of turbulent 
diffusivity (more precisely, the coefficient of the second 
order spatial derivative term) is not flijk but is instead 



(10) 



As a complementary approach, which reinforces the 
above arguments, we consider, via a classical perturba- 
tion analysis, the evolution, under a small-scale velocity 
field, of a magnetic field with a long-wavelength modula- 
tion; this is the approach adopted by Roberts Q in his 
analysis of two-dimensional cellular flows. We write the 
mean field in the form Bo = Bq exp(ik ■ x + pi) and the 
fluctuating field as b = bexp(ik ■ x +pt), where b varies 
on the same small spatial and temporal scales as u. The 
wavenumber k = \k\ is assumed to be small; accordingly, 
we may develop the solution in powers of k] specifically, 
the growth rate has the expansion 



p = kpi + k 2 p 2 



(11) 



The coefficient p 2 gives an unambiguous description of 
the diffusion term. 

At leading order, the fluctuating and mean induction 
equations take the form 

{0 t - V V 2 )b = B ■ Vm + V x G , (12) 
piB = iky. £ Q ; (13) 

and at the next order, 

(d t - r/V 2 )bi + p 1 b Q - 2ink ■ Vb = 

-i(u ■ k)B + V x Gi + ik x G , (14) 
p 2 B = ikx£i - r)B (15) 

(where we have assumed V • u = for simplicity). The 
formal means of solution is now clear, although the in- 
version of the various linear operators may of course lead 
to technical difficulties. The fluctuating field bo can be 
obtained from (TT2"]) ; £o can then be evaluated and hence, 
from (fT3j) . pi determined as 



Pi 



ik x £q ) B 



B 2 



(16) 



Proceeding, equation (TT4"]) can be solved (formally) for 
b%; £\ can then be evaluated and hence p 2 calculated, 
from (fl~5|) . as 



P2 = 



(ik x£i- r/B ^ ■ B* 



(17) 



The traditional diffusive term (i.e. that represented by 
the (3 tensor) comes not from solving (TT4")) . but instead 
from solving the equation 

(d t -ryV 2 )Si = ~i(u-k)B + V x G^ (18) 
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It is clear that this does not tell the whole story regarding 
the growth rate; the term pib in (fT4")l . the effects of 
which are encapsulated in the T tensor, may be equally 
important. 

We may be more explicit if we consider a simple ex- 
ample in which u (and hence b) takes the form u = 
M.(u exp(i(m ■ x — u>t)) . The awkward G terms then van- 
ish, allowing explicit solution of equations (TT2"|) and (fTi)) 
for all Rm. Then £q is given by 



e 



X R 1 — X (' m '-^o)" 
2 I V 



(19) 



where D = —iuj + r/m 2 . The quantity pi then follows 
from expression (|16p . At the next order we obtain 



Si 



_1 | m* x (ft ■ fc)J3 
2 1 V 

2rj(m ■ k)u* x (im ■ Bq)u 
V 2 

Piu* x (im • JBo)w 1 

[ ' 



(20) 



with p2 given by expression (|17[) . The first two terms 
are encapsulated in the traditional fiijk term. However, 
the final term, being proportional to pi, can be identi- 
fied with the first term in (|10[) . This term arises from 
the fact that we are solving an eigenvalue problem for 
the growth rate, rather than studying the statistically 
steady response to an imposed current. We note here 
that the final term in expression (|20[) is small compared 
with the other terms when Rm is small. If for example 
we regard r\ as fixed and thus take Rm oc \u\, we see that 
the last term oc p\Rm 2 ~ Rm 4 , while the other terms 
oc Rm 2 when Rm is small. In this situation the First 
Order Smoothing Approximation can be used for quite 
general velocity fields, leading to formulae analogous to 
(|19l20p . However when Rm is not small, the new term 
may be at least as large as the others, while direct cal- 
culation is not possible. 

It is therefore of interest to investigate, for a specific 
flow, the relative magnitudes of the two contributions 
to the diffusivity in expression (|10[) . particularly when 
Rm is much greater than unity. This can be most read- 
ily achieved by considering spatially-periodic velocities 
of the form u(x,y, i); for such flows, the magnetic field 
takes the form B(x, t) = B(x, y) exp(pt + ikz). By con- 
sideration of the kinematic dynamo problem for a range 
of small values of k, one can obtain the first two terms 
in the expansion for the growth rate, 



p = p\k +P2k 2 



(21) 



The real part of p\ provides information on a^ , the real 
part of P2 encompasses all of the diffusive contributions. 
Via the complementary approach of imposing two inde- 
pendent spatially uniform magnetic fields of the form, for 



example, 

B = {B + C t, 0, 0) and B Q = (0, B + C t, 0), (22) 

constraining the perturbation magnetic field to be inde- 
pendent of z (thus eliminating the possibility of small- 
scale dynamo action), and then taking spatial averages 
over the iry-plane and temporal averages over the fast 
flow time scale, we may independently calculate the ten- 
sors otij and Tij (here 1 < i,j < 2). We are then in a 
position to quantify the two contributions to the mag- 
netic diffusivity in expression (fT0|) . We should note that, 
at least in theory, it is also possible to obtain an inde- 
pendent measure of the elements of the flijk tensor, by 
calculating the emf after the imposition of steady mag- 
netic fields with a uniform gradient. However, this is ex- 
tremely computationally expensive, in comparison with 
the calculations of the a and T terms, since it requires a 
spatial domain that is very large compared with a typi- 
cal velocity scale; by contrast, the computational domain 
needed to calculate the a and T terms is simply that of 
the spatial periodicity of the flow. 

Here we consider the specific two-dimensional time- 
dependent flow introduced by Otani [5fl (sometimes re- 
ferred to as the MW+ flow), given by 

u(x, y,t) = V x (ip(x, y, t)z ) + ip(x, y, t)z , (23) 

where 



ip(x, y, t) — 2 cos 2 t cos x — 2 sin 2 1 cos y. 



(24) 



Symmetry considerations @, 0] dictate that onj and Pijk 
take the form = a&y, Pijk = fieijk; similarly it can 
be shown that = TSij . Thus we need only consider 
an imposed field in either the x or y-direction in order 
to determine ctij and r^ . It follows from (TlU)) that the 
'new' contribution to the diffusivity is — aT. 

Figure DJa) shows the spatially, but not temporally, 
averaged emi £ x (i.e. (uxb) x ) versus time, for an imposed 
field B = tx and for Rm = 100. Figure [Tfb) shows B x , 
the emf after a further temporal average over the fast 
time scale of the flow. Since, for an imposed field of the 
form Bo = (Bo + Cot)x, £ x takes the form 



£ x =a(Bo + C t) + rC'o, 



(25) 



a and T can be evaluated from the slope and intercept of 
the straight line in Fig. [T{b). (Note that we can obtain 
both quantities simply by having Co non-zero; setting 
Bo 7^ 0, Co = (the traditional means of computing a) 
supplies an independent verification of a). 

Table U shows, for a range of Rm, the values of p\ and 
P2, obtained from solving the kinematic dynamo problem 
for a range of small values of k, together with the values 
of a and T, determined from the emf after the imposition 
of a spatially uniform field Bo = tx. The final column 
contains the —aT contribution to the diffusivity; the full 
diffusivity is given by — p2 = — aT. (By the very nature 
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time 

FIG. 1: (a) Spatially averaged emf, £ x , versus time for the 
MW+ flow, for Rm = 100 and an imposed field Bo = tx. 
(b) £ x , obtained from over the fast time scale of 

the flow, versus time. 

of the means of determining these values, with the pj ob- 
tained from fitting a parabola to the p versus k curve, 
whereas a and T are determined from long temporal av- 
erages, the former (particularly^) cannot be determined 
to the same level of accuracy as the latter.) 

An interesting pattern emerges from Table HI At low 
values of Rm the diffusivity is predominantly due to [3; 
this is consistent with the explicit result for the simple 
monochromatic flow discussed above. However, as Rm is 
increased, the —aT term increases in importance until, at 
Rm = 100 it is the dominant contribution. For Rm = 50 
and 100 the value of (3 is negative, though small in magni- 
tude in comparison with the total diffusivity — pi- This is 



certainly possible, though it should also be borne in mind 
that the determination of pi involves some slight inaccu- 
racy. The crucial element though, which is unaffected 
by small changes in p2 , is that the magnetic diffusion at 
these higher values of Rm comes almost completely from 
the — aT contribution. 

To summarize, we have identified an inconsistency in 
the usual means of expressing the mean emf only in 
terms of spatial derivatives. In the general expression 

TABLE I: pi and p2, calculated from the p vs k curve after 
solving the kinematic dynamo problem, together with a and 
r, calculated from consideration of the emf after the impo- 
sition of the mean field Bo = tx, for a range of values of 
Rm. 



Rm 


Pi 


P2 


a 


r 


- aT 


1 


0.80 


-1.56 


-0.80 


0.20 


0.16 


10 


0.75 


-0.93 


-0.75 


0.85 


0.64 


50 


1.51 


-6.2 


-1.52 


4.4 


6.69 


100 


1.9 


-14.5 


-1.92 


7.8 


14.98 



for the emf it is vital to retain temporal derivatives of 
the mean field; substitution in terms of spatial deriva- 
tives can be performed, but this needs to be done in a 
consistent order-by-order manner. At small values of Rm 
the new term contribution to the magnetic diffusivity is 
small; however, in the example we have studied, it be- 
comes dominant at large Rm. 

For clarity we have focused in this paper on a sim- 
ple mean field problem with no mean flow. We could 
straightforwardly have included this term in the equa- 
tions, which would have changed the dynamo growth rate 
and hence the value of the turbulent diffusivity term (p2). 
A consequence of the work described here is that even if 
the mean flow can be ignored in the fluctuating equa- 
tions, it may still affect the diffusivity since this is linked 
to the growth rate. 
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